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Abstract. Two-armed, grand design spirals and inner and outer rings in barred galaxies 
can be due to orbits guided by the manifolds emanating from the vicinity of the L\ and L 2 
Lagrangian points, located at the ends of the bar. We first summarise the necessary theoreti- 
cal background and in particular we describe the dynamics around the unstable equilibrium 
points in barred galaxy models, and the corresponding homoclinic and heteroclinic orbits. 
We then discuss two specific morphologies and the circulation of material within the corre- 
sponding manifolds. We also discuss the case where mass concentrations at the end of the 
bar can stabilise the L t and L 2 and the relevance of this work to the gas concentrations in 
spirals and rings. 
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1. Introduction 

Barred galaxies often have interesting mor- 
phological characteristics. These include two- 
armed, grand design spirals and rings, both 
inner and outer. Over the last five years we 
developed a theory which accounts for the 
formation and the properties of these struc- 
tures, starting with the dynamics of the re- 
gions around the ends of the bar. Our main re- 
sults are summarised in four papers [Romero- 
Gomez et al. 2006 (Paper I), 2007 (Paper 
II), 2008; Athanassoula, Romero-Gomez & 
Masdemont 2008 (Paper III)], while a fifth one 
focuses more on comparisons with observa- 
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tions (Athanassoula, Romero-Gomez, Bosma 
& Masdemont 2009, Paper IV). In these pa- 
pers, the reader will also find other relevant ref- 
erences, concerning other aspects of this prob- 
lem, either theoretical, or observational. 

As discussed in Binney & Tremaine 
(2008), barred galaxies have five equilibrium 
points, often referred to as Lagrangian points, 
of which two are unstable {L\ and L2) and 
three are stable (L3, L4 and L5). The former are 
located on the direction of the bar major axis, 
while the latter are located at the centre of the 
galaxy and on the direction of the bar minor 
axis. From the region around the L\ and L2 
emanate the manifolds which guide the chaotic 
orbits escaping this vicinity. These manifolds 
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can be simply thought of as tubes that guide 
and confine the chaotic orbits in question. It 
is this confined chaos that can account for the 
grand design spiral arms and the inner and 
outer rings. 

The structures we wish to understand are 
the spirals and the inner and outer rings 
in barred galaxies. Inner rings have roughly 
the size of the bar and are elongated along 
it. They are noted as r in the classifica- 
tion of de Vaucouleurs and his c ollaborators 
dButa. Corwin & Odewahn 1 120071 and refer- 
ences therein). Outer rings are similar, but of 
larger size. They are either elongated perpen- 
dicular to the bar (R\ type), or along it (^2 
type). In some cases the two types of outer 
rings co-exist, and the structure is then noted 
R\R2- Thus, a barred galaxy with both an inner 
ring and an outer ring of the first type will be 
noted rR\. A description of this classification, 
as well as more inform ation on the se struc- 
tures has been given by Buta d 19951) . Spirals 
are also very common amongst barred galaxies 
dElmegreen & Elmegreenll989l) . In most cases 
they start off from the ends of the bar, and they 
always wind outwards in a trailing sense. 

This paper is organised as follows. In 
Sect. |2] we describe the bar models used in our 
calculations. In Sect. [3] we present the theoret- 
ical basis of our theory. In Sect. [4] we discuss 
two specific morphologies, the spirals and the 
rR\ rings, and the different patterns of matter 
circulation that they entail. We also discuss the 
types of potentials that can lead to these mor- 
phologies. In Sect. [5] and [6] we address two 
physical problems, namely the presence of two 
extra mass concentrations centred on the two 
ends of the bar and the effect of the gas, re- 
spectively. Finally, we conclude in Sect. [7] 

Response calculations in barred galaxy po- 
tentials were discussed in this meeting by 
Patsis i n his t alk (these proceeedings; see also 
iPatsisI 42006) and references in both), while 
Efthymiopoulos et al. put up a poster dis- 
cussing dri ving by spirals (these proceedings; 
see a lso iTsoutsis. Efthymiopoulos & Voglisl 
(2008) and references in both). 



2. Models 

We model the barred galaxy as the superposi- 
tion of an axisymmetric component with a non- 
axisymmetric one, representing a bar rotating 
around the centre of the galaxy. 

We use three different bar models to de- 
scribe the bar component. The basic bar model 
is de scribed by a Ferrers ellipsoid dFerrersI 
1 1 877b of density distribution 



where m 2 = x 2 /a 2 + y 2 lb 2 . The values of a 
and b determine the shape of the bar, a being 
the length of the semi-major axis and b being 
the length of the semi-minor axis. In a frame 
of reference co-rotating with the bar, the major 
axis is placed along the x coordinate axis. The 
parameter n measures the degree of concentra- 
tion of the bar, while the parameter po repre- 
sents the central density of the bar. For these 
models, the quadrupole moment of the bar is 
given by the expression 

Q m = M b {a 2 - b 2 )/(5 + 2b), 

where Mb is the mass of the bar, equal to 

M b = 2 {2n+3) nab 2 p r(n + l)F(n + 2)/T(2n + 4) 

and F is the gamma function. 

This bar model is the basic model of 
lAth anassou ial d!992al) and it will also be our 
basic bar model. We refer to it as model A and 
it has essentially four free parameters that de- 
termine the dynamics in the bar region. The ax- 
ial ratio a /b and the quadrupole moment (or 
mass) of the bar Q m (or Mb) will determine 
the strength of the bar. The third parameter is 
the pattern speed, i.e. the angular velocity of 
the bar. The last free parameter is the central 
density of the model (p e ). For reasons of con- 
tinuity we will use the same nu merical values 
for the model parameters as in lAthanassoulal 
dl992allbh and in Papers II and III. The length 
of the bar is fixed to 5 kpc. Mo re information 
on thes e models can be found in Athan assoula! 
dl992al) and in Papers II and III. 

Since we want to study the formation of 
spiral arms and rings and both features are lo- 
cated in the outer regions of the bar or well 
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outside it, we consider two other bar poten- 
tials to avoid a well known disadvantage of 
Ferrers potential, namely the fact that the non- 
axisymmetric component of the force drops 
very steeply beyond a certain radius, particu- 
larly for values of the density index n larger 
than one (eq.[T]i, so that the axisymmetric com- 
ponent dominates in the outer regions. These 
ad-hoc potentials are of the form eA(r) cos(20), 
e being related to the strength of the bar. 

The first ad-hoc bar po tential we use is 
adapted from lDehnenl J2000) and has the form 



$(r,#) = --ev 2 cos(26») 



2 - (r/af , r<a 



(a/r) n , 



r > a. 



The parameter a is a characteristic length 
scale of the bar potential and Vrj is a constant 
with units of velocity. The parameter e is a free 
parameter related to the bar strength. In this pa- 
per we use a — 5 and n - 0.75. We will refer 
to models with this non-axisymmetric compo- 
nent as the D models. 

Finally, our third model has the bar poten- 
tial: 

0>(r, &) = e Vr(ri - r) cos(26>), 

where r\ is a characteristic scale length 
of the bar potential, which we will take 
for the present purposes to be equal to 
20 kpc. The parameter e is related to the 
bar strength. This type of model has al- 
ready been widely used in studies of bar 
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dynamics (e.g. iBarbanis & Woltierl 
Contopoulos & Papayannopoulosl 
Contopoulosl 1 1 98 lb . We will refer to models 
with this non-axisymmetric component as the 
BW models. 

Models D and BW are meant to represent 
and include forcings not only from bars, but 
also from spirals, from oval discs and from tri- 
axial haloes. We have kept all these forcings 
bar-like, i.e. we have not included any radial 
variation of the azimuthal dependence. This 
was done on purpose, in order not to bias the 
response towards spirals and in order to avoid 



the extra degree of freedom resulting from the 
azimuthal winding of the force. 

Throughout this paper we use the follow- 
ing system of units: For the mass unit we take 
a value of 1O 6 M , for the length unit a value 
of 1 kpc and for the velocity unit a value of 
1 km/sec. Using these values, the unit of the 
Jacobi constant will be 1 km 2 /sec 2 . 

3. Theoretical background 

3.1. Equilibrium points and dynamics 
around them 

We will now study the dynamics of the motion 
in the equatorial plane of the galaxy, z — 0. We 
work in a frame of reference co-rotating with 
the bar, i.e. a frame in which the bar is at rest 
and by convention we place it along the x axis. 

The energy of a particle in this rotating 
frame is a constant of the motion and is given 
by 

Ej = -{x 2 + f) + ^ff, 

where <£ eff = O - ±£2 2 (x 2 + y 2 ) is 
the effective potential, i.e. the potential in 
the rotating frame, and O is the potential 
dBinney & Tremainell2008h . The curve defined 
by ^eff = Ej is the zero velocity curve (ZVC) 
and divides the equatorial plane in three differ- 
ent regions, namely the forbidden region, the 
inner region, where the bar is located, and the 
outer region (see Fig. [TJ. The forbidden region 
is enclosed within the ZVC and is forbidden to 
particles with energy equal to that of the ZVC, 
or smaller. The inner region is delineated by 
the inner part of the ZVC of energy equal to 
that of the L\ and L2 and includes the centre, 
while the outer region is delineated by the outer 
part of this ZVC and includes the outermost 
parts of the galaxy. 

The system has five equilibrium points lo- 
cated on the equatorial plane where 



SO, 



eff 



d0>, 



elf 



dy 



= 0. 



L\ and L2 lie on the direction of the bar major 
axis and they are saddle points, i.e. linearly un- 
stable. The other three are stable. L3 is at the 
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Fig. 1. Location of the equilibrium points (L,, i = 
1, .., 5, marked with asterisks) in model A with 
a/b=5, Q m = 4.5 x 10 4 , r L = 6 and p c = 2.4 x 10 4 . 
The dashed line gives the outline of the bar, while 
the dot-dashed lines show the zero velocity curves 
for a characteristic energy. Around L t and L 2 we plot 
the Lyapunov periodic orbits of this energy (thick 
solid lines). 

origin of coordinates and it is a minimum of 
the effective potential, while L\ and L5, located 
along the direction of the bar minor axis (see 
Fig. [TJ, are maxima. We define as the dis- 
tance from the origin to L\(L 2 ). 

Around each of the equilibrium points, 
there exists a family of periodic or- 
bits. For example, around L3 we have 
the X] family of stable perio dic orbit s 
dConto noulos & Pa pavannopoulosl [l980), 
which is k nown to be responsible fo r the bar 
structure (Athanassoula etaD Il983l) . In the 
linear approximation, the motion around L\ 
and L2 consists of the superposition of an 
exponential part and an oscillation: 

x(f) = Xie At + X 2 e~ M + X-i cos(wf + </>), 
y(t) = A x X x e M - A x X 2 e~ M 
+ A 2 X } sin(wf + <p), 

where +A, ±a>i are the eigenvalues of the dif- 
ferential matrix at L\(L2), and Xi,X 2 ,X-i,A\, 
and A 2 are constants. For each energy level 
there exists a unique periodic orbit around L\ , 
or L 2 , called Lyapunov orbit (lLyapunovll l949). 



To obtain them, we take initial conditions from 
Eq.OwithXi = X 2 = 0: 



x(f) = (x,y,x,y) = 

(X3 cos(w/ + <p), 
A2Z3 sin(w? + <p), 
-X^tosmiut + <p), 
A2X3C1J cos(u)t + <p)). 



These orbits can of course be calculated 
numerically directly from the equations of mo- 
tion in the adopted potential and reference 
frame, without the use of the linear approxi- 
mation. Examples are shown in Fig. Q] They 
are unstable, so they cannot trap particles like 
a stable family. However, there are other ob- 
jects emanating from the periodic orbits that 
can trap particles. They are called the invari- 
ant manifolds and are associated to the periodic 
orbits. From each periodic orbit emanate four 
branches of asymptotic orbits, two of them be- 
ing stable and the other two unstable. The char- 
acterisation of stable and unstable in this def- 
inition is in the sense that the stable invariant 
manifold is a set of orbits that approach asymp- 
totically the periodic orbit, while the unstable 
invariant manifold is a set of asymptotic orbits 
that depart asymptotically from the periodic or- 
bit. That is, there are four favourite directions 
along which material can approach or escape 
from the vicinity of the corotation region and 
they are given by the stable and the unstable 
manifolds, respectively. 

In the linear case, we can obtain an initial 
condition for the stable invariant manifold con- 
sidering X\ — and X 2 + in Eq. 13.11 The 
exponential term proportional to X 2 vanishes 
when time tends to infinity and the trajectory 
tends to the periodic orbit. Analogously, if we 
consider initial conditions with X\ + and 
X 2 - 0, we will obtain the unstable invariant 
manifold. 

Thus, the manifolds can be thought of as 
tubes which guide chaotic orbits escaping the 
vicinity of L\ and L 2 . Therefore, according to 
this theory, spiral arms will be due to confined 
chaos. 
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3.2. Transfer of matter. Homoclinic and 
heteroclinic orbits 

Invariant manifolds are not restricted to the 
neighbourhood of L\(L,2), but they extend 
well beyond the region around the equilibrium 
points. Two of the branches, one stable and one 
unstable, are located in the inner region, while 
the other two are located in the outer region, 
so that the invariant manifolds connect the two 
regions and the periodic orbits act as a gate- 
way. In Fig. 12 we can see four examples of 
the morphologies obtained with these invariant 
manifolds. 

The stable (respectively, unstable) mani- 
folds of a periodic orbit cannot intersect in the 
phase space. However, the stable and the unsta- 
ble branches can intersect each other. If these 
stable and unstable branches are associated to 
two different periodic orbits, one around the L\ 
and the other around the L2, we obtain hetero- 
clinic orbits, i.e., orbits that connect one of the 
ends of the bar with the opposite one. Due to 
the symmetry of the system, the particles fol- 
lowing these orbits will outline a morphology 
similar to that of an R\ ring (left panel of Fig. 

If, on the other hand, the stable and un- 
stable branches are associated to the same 
Lyapunov periodic orbit, the intersection gives 
rise to what is known as homoclinic orbits. 
That is, asymptotic orbits that connect one of 
the bar ends with itself. Considering the sym- 
metry of the system, we will have particles out- 
lining trajectories reminiscent to that of R\Ri 
rings, i.e. the particles will form two outer 
rings, one with major axis perpendicular to the 
bar major axis {R\ ring) and one with major 
axis parallel to it (R2 ring), as shown in the sec- 
ond panel of Fig. [2] 

If the stable and unstable manifolds ema- 
nating from one of the Lyapunov orbits do not 
return either to that orbit, or to the correspond- 
ing one around the Lagrangian point at the op- 
posite side of the bar, but unwind outwards, 
they form a spiral shape (third panel of Fig. 

Thus, such morphologies can give rise to 
two trailing grand design spirals. We will call 
orbits following such manifolds escaping, be- 



cause they can escape the vicinity of the bar 
and reach the outer parts of the galaxy. 

Finally, if the outer branches of the unsta- 
ble invariant manifolds emanate with an ap- 
propriate pitch angle, they will intersect in 
the configuration space forming an outer ring 
whose major axis is parallel to the bar major 
axis, i.e. an R2 ring (right panel of Fig. [2). 

4. Two specific examples 

As described in the previous section, manifolds 
can reproduce the right morphologies for both 
spirals and rings (inner, R\, R2, RiRi)- In this 
section we will consider in more detail two 
possible morphologies, rRi rings and spirals. 

An rRi morphology necessitates hetero- 
clinic orbits. In this case, the stable and un- 
stable manifolds overlap in the configuration 
space. The inner branches form an inner ring, 
elongated along the bar major axis. The outer 
branches form an 8-like or 0-like shape, i.e. an 
R\ ring morphology. In particular, this is elon- 
gated perpendicular to the bar major axis (i.e. 
its major axis is perpendicular to that of the 
bar) and it also produces characteristic dimples 
on either side of the bar, at or near the ends 
of the inner ring. Such dimples have also been 
fou nd, in the same locatio n, in observed galax- 
ies (lButa&Crockerlll99ll) . 

This type of manifold morphology leads 
to interesting circulation patterns. Particles can 
follow essentially four different paths, shown 
in Fig. [3] One possibility is that they circu- 
late along the inner branches of the invariant 
manifolds forming the inner ring (left panel 
of Fig. [3). This circulation is anti-clockwise. 
The second possibility is that they follow the 
outer branches of the manifolds, outlining the 
outer ring (middle panel of Fig. |3). This cir- 
culation is clock-wise. Finally, they can have 
a mixed trajectory, i.e. follow both inner and 
outer branches. Thus, particles starting from 
L\ and following the inner branch of the un- 
stable manifold approach the neighbourhood 
of L2, move out of the bar region and follow 
the outer branch of the unstable manifold. They 
then reach a maximum distance from the cen- 
tre and come back to the neighbourhood of L\ 
(right panel of Fig. |3). In this way the circu- 
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Fig. 2. Four different morphologies. First panel: Model with heteroclinic orbits having an rRi morphology. 
Second panel: Model with homoclinic orbits having an R]Ri ring morphology. Third panel: Model with 
no heteroclinic or homoclinic orbits having two spiral arms. Fourth panel: Model with no heteroclinic or 
homoclinic orbits having an R2 ring. 
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Fig. 3. The four possible circulation patterns for particles on orbits guided by the manifolds. Any given 
such particle can stay on the inner branches of the manifold, i.e. can trace the inner ring (left panel), or it 
can stay on the outer branches of the manifold, i.e. can trace the outer ring (middle panel). Alternatively, it 
can trace a path that includes one inner and one outer branch. Two such paths are possible. One is shown 
in the right panel. The other, i.e. the fourth alternative, is not plotted here since it is similar to the third one, 
except that it is symmetric with respect to the bar major axis. After a particle has completed a full circulation 
path, it can either repeat the same path, or take any of the other three ones. The manifolds coming into play 
in each case are plotted in black in the corresponding panel, and the remaining ones in light grey. In white, 
we plot the Lyapunov periodic orbits of the corresponding energy level. The arrows show the direction of 
the motion. 



lation pattern is completed. A fourth circula- 
tion pattern (not shown in Fig. [3} has an iden- 
tical shape but is reflected with respect to the 
bar major axis. In both the third and the fourth 
cases the circulation is clockwise with respect 
to their 'centres', i.e. with respect to the L3 and 
L4 Lagrangian points, respectively. 

There are, thus, four circulation paths in 
total, one using only the inner branches, one 
using only the outer branches and two mixed 
ones, using both inner and outer branches. 



Once one such pattern is completed, a particle 
can either repeat it or take any of the other three 
circulation patterns. Thus, material from inside 
corotation can move outside it and vice-versa. 
However, the maximum radius such material 
can reach is bound by the maximum distance 
of the path from the centre, located on the di- 
rection of the bar minor axis. 

As a second example, let us discuss a spi- 
ral morphology. Material, initially on an in- 
ner stable manifold branch, can, via the neigh- 



Athanassoula et al.: Rings, spirals and manifolds 



bourhood of one of the unstable Lagrangian 
points, move to an outer unstable branch. If this 
is of the escaping type (see previous section), 
we obtain a morphology similar to that of the 
grand design spiral arms of a barred galaxy. 

The circulation pattern in this case is much 
simpler (Fig. 0J. Matter from the inner region 
(more specifically from the outer regions of the 
bar or its immediate vicinity) moves along the 
inner manifolds towards one of the two un- 
stable Lagrangian points and from there onto 
the corresponding outer branch of the unstable 
manifold. It can thus escape the inner region 
and reach outer parts of the disc. 

Of course, an inwards going route could 
also be possible, at least in principle. Then mat- 
ter from the outer parts of the disc could move 
inwards along a stable outer manifold branch 
to the L\ or the L2 and from there into the in- 
ner region. However, as will be discussed fur- 
ther in Sect. [7] the existence of a given man- 
ifold does not necessarily imply the existence 
of the corresponding structure in a galaxy, or in 
an A^-body simulation. For this, the manifolds 
have to confine a sufficient amount of matter. 
This is similar to periodic orbits, where the ex- 
istence of a stable family will not imply the for- 
mation of any structure if it does not trap any 
regular orbits around it. Thus the route bring- 
ing matter inwards, although in principle possi- 
ble, would imply the existence of leading two- 
armed grand design spirals in barred potentials 
and has not yet been observed either in simula- 
tions or in real galaxies. Whether and how mat- 
ter gets trapped by periodic orbits, or by man- 
ifolds depends on the formation history of the 
galaxy and is well beyond the scope of this pa- 
per. In Paper IV, however, we examine a few 
specific cases. 

The two different morphological types dis- 
cussed above, namely rR\ rings and two-armed 
grand design spirals, do not come in the same 
potentials. The rR\ morphology comes in mod- 
els with bar potentials which are not too strong 
in the corotation region and immediately be- 
yond it. In Paper III we gave, for the types 
of models of Sect. [2] upper limits of the bar 
strength beyond which this morphology will 
not occur. On the other hand, spirals form in 
models with stronger bar or appropriate spiral 
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Fig. 4. Circulation pattern for particles in models 
with spiral arms. The particle will move along the 
inner stable manifold branch and reaching the vicin- 
ity of one of the two unstable Lagrangian points will 
continue outwards along the unstable branch of the 
manifold. The direction of the motion is given by the 
arrows. We also plot the ZVC of the same energy as 
the manifolds (solid line) and give the position of 
the Lagrangian points with asterisks. 



forcings. We also showed that it is the strength 
of the non-axisymmetric potential at and im- 
mediately beyond corotation that are the best 
indicators of the morphological type. For this 
we use the quantity 



Q,(r) = (dO(r, 0)/d0) max /(rd®o/dr), (2) 



where <L> is the potential, <E>o is its axisymmet- 
ric part and the maximum in the numerator is 
calculated over all values of the azimuthal an- 
gle 8. We calculate Q t (f) at the radius of L\, 
i.e. at r — rn and denote it by g/,Li ■ As shown 
in Paper III, R\ morphologies form in models 
with relatively low values of Q, ^ , while higher 
values of this quantity give rise to other mor- 
phologies. 
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5. Cases with stable L\ and L 2 
Lagrangian points 

5.1. Achieving stability 

The morphologies and the characteristics of 
the invariant manifolds described so far as- 
sume that the L\ and L2 equilibrium points are 
unstable and with linear stability of the type 
"saddlexcentrexcentre". We will hereafter re- 
fer to this case as the standard case. To make 
our theory more general, we also examine un- 
der what conditions the L\ and L2 can become 
stable and what consequences this will have on 
the galaxy morphology (Paper III). 

It is possible to stabilise the equilibrium 
points by adding a concentration of mat- 
ter around them. We tested this by adding 
two identical, small Kuz'mi n/Toomre discs 
dKuz' mini 1956t iToomrdl 1963b , one centred on 
each of the L\ and L2. By increasing their mass, 
or decreasing their scale-length, the equilib- 
rium point L\ bifurcates, becoming stable, 
while two new unstable points appear, both on 
the direction of the bar major axis. They are lo- 
cated one on either side of L\, and are called L' x 
and L", the subscript, i or o, denoting whether 
the corresponding new Lagrangian point is lo- 
cated inside or outside the L\. Analogously, 
L2 also becomes stable and two more unstable 
equilibrium points, L' 2 and L° y appear on either 
side of it. The mass of the added concentration 
at which this is attained is a decreasing func- 
tion of its scale-length. Similarly, the concen- 
tration scale-length at which this is attained is 
an increasing function of its mass. 

This new configuration is shown in Fig. [5] 
where we plot the ZVC for a given energy, the 
position of the equilibrium points and the out- 
line of the bar, for model A with ajb = 2.5, 
r L = 6, Q m = 4.5 x 10 4 , and p c = 2.4 x 
10 4 . The scale-length and the mass of the two 
identical Kuz'min/Toomre discs are 0.6 and 
0.025M/,, respectively. The linear stability of 
the four new equilibrium points is of the type 
"saddlexcentrexcentre", so we can compute 
the family of unstable periodic orbits and the 
manifolds associated to them, as we did for 
those of L\ and Li in the standard case. 
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Fig. 5. Location of the equilibrium points (marked 
with an asterisk) in an A model with a/b = 2.5, 
n = 6, Q m = 4.5 x 10 4 and p c = 2.4 x 10 4 , to which 
have been added two identical Kuz'min/Toomre 
discs centred at the L\ and L2 Lagrangian points. 
The scale length and the mass of these discs are 0.6 
and 0.025Mj, respectively. The dashed line marks 
the position of the bar, while the dot-dashed lines 
show the zero velocity curves. 



5.2. Global morphology 

In Fig.|6]we compare the morphology obtained 
with the standard case to that with the same 
bar model plus the two mass concentrations 
around the L\ and L2 in the form of two small 
Kuz'min/Toomre discs. The two morphologies 
are essentially of the same type, since in both 
cases we obtain an rR\ structure. However, the 
sizes and axial ratios of the two rings change 
drastically when the two mass concentrations 
are added and the number of equilibrium points 
increases to 9. The inner ring becomes smaller 
and is located very near the outline of the bar. 
Its axial ratio approaches that of the bar, while 
in the 'standard' case it has a smaller elliptic- 
ity. The major axis of the outer ring becomes 
much larger, and its axial ratio changes accord- 
ingly. If no material circulates between L\ and 
L° and between L' 2 and L° v the inner and outer 
rings will not join, i.e. the outer ring will be de- 
tached. Such detached rings are often observed 
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Fig. 6. Comparison of the manifold morphology in the standard model with 5 equilibrium points (left panel, 
equilibrium points as in Fig.lTJ and with 9 equilibrium points (right panel, equilibrium points as in Fig. [5). 
The bar outline is given by a dashed line and the corresponding Lyapunov orbits by a black solid line. 



and this may be the way, or one of the ways, to 
explain their formation. 

5.3. Specific morphology around the L\ 
and L 2 

The morphology of the unstable and stable 
manifolds in the immediate neighbourhood of 
the L\ Lagrangian point is shown in Figs.|7]and 
[8] respectively, for a model which has a spiral 
morphology. This has several similarities, but 
also several differences from the model shown 
in Figs. 14 and 15 of Paper III, which is for a 
less strong bar with manifolds of rR\ morphol- 
ogy- 

The Lyapunov periodic orbits around L\ 
are elongated in the same direction as the bar, 
while the periodic orbits around L\ and L" 
are elongated perpendicular to it. It is the one 
around L\ that is the largest and the one around 
Lj that is the smallest, the third one being of in- 
termediate size. 

There are now in total four unstable 
branches (Fig. |7}. The outer branch of the L° 
manifold, as well as the inner branch of the 
Lj have a very simple morphology, while the 
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Fig. 7. Morphology of the unstable manifold 
branches in the vicinity of a stable L\ Lagrangian 
point (see text). The two upper panels show the man- 
ifolds linked to L". and the two lower ones those 
linked to L\. The arrows show the direction of the 
motion around the manifolds and the dashed line 
shows the outline of the bar. Three periodic orbits 
- one around each of the L\, L\ and L" points - are 
given with solid lines. 
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Fig. 8. Same as Fig. [7] but for the stable branches 
of the manifolds. 



two other branches have a more complicated 
and interesting morphology. The inner branch 
of the L° emanates from the vicinity of L° and 
circumventing the L\ from the above extends 
towards L\, but before reaching it turns back 
towards L ., now circumventing the L\ from 
below. The unstable outer branch of L\ goes 
towards L" circumventing L\ from below and 
then moves outwards beyond the corotation re- 
gion. The morphology of the stable branches is 
similar and is shown in Fig. [8] 



6. Ersatz for gas 

Observations show that gas is intimately 
linked to both spirals and rings. For exam- 
ple, galaxies with smooth, very gas poor 
spiral arms, known as anemic, are a rarit y 
rivandenBerghlll976t lElmegreen et al.ll2002l) . 
Furthermore, the formation of rings and 
spirals from gas has been witnessed in a 
(e.g. Schwar3 



number 



of simulations 
n 



1981 



Combes & Gerin 1985; Athanassoula 1992b; 
Lindblad, Lindblad & Athanassoula 1996; 
Wada & Kodal 1200 it iLin. Yuan & Butal mM. 

~W). 



Rautiainen, Salo & Laurikainen 2008). We 



should thus compare the dynamics of the gas 
with that of the manifolds presented here. 
Gas, however, has different equations of 



motion from stars, so we need to modify our 
calculations accordingly. 

ISchwarzl (I1979L 1 19811 11984 1 19851) uses 
sticky particles to simulate the gas and mod- 
els collisions in a particularly straightforward 
way, so we can introduce a similar procedure 
also in our calculations. In Schwarz's simula- 
tions, particles represent gaseous clouds which 
lose a certain fraction of their kinetic energy 
when they collide. In practise, the y lose a 
certai n fraction / of their velocity dSchwara 
1981), or only of the component of this ve- 
locity tha t is along the l ine j oining the two 
particles dSchwarzll 19841 Il985l) . Thus, v 2 = 
±(1 - f)\'\, where + is for the tranverse veloc- 
ity component, - is for the component along 
the line joining the two particles, and v\ and v 2 
are the velocities before and after the collision, 
respectively. The values off in Schwarz's sim- 
ulations range between 0.8 and 1. 

We introduce collisions in our simulations 
in a similar way, thus obtaining an ersatz for 
gas, as described in more detail in Papers 
III and IV. We also use the same models as 
Schwarz in order to be able to make compar- 
isons. Here we present results for his standard 
model with a = 2.6, q = 0.1, £l p = 0.27 and/ 
0.8. for which he gives sufficient results and in- 
formation in his papers to allow comparisons, 
and which has an rR\ manifold morphology. 
In this model we calculated a number of or- 
bits of the outer branch of the unstable mani- 
fold and drew random numbers to find the lo- 
cations along its trajectory where each particle 
will undergo a collision. In the particular ex- 
ample illustrated here we have three collisions 
per half bar rotation. We tried, however, dif- 
ferent numbers of collisions, around the values 
used by Schwarz, and found qualitatively sim- 
ilar results. 

To determine the result of a collision, we 
take a small box around the collision position 
and calculate the average velocity of all or- 
bits in that box. We then decrease the veloc- 
ity of the particle relative to that mean by a 
factor +(1 - /). Fig. |9]shows a comparison of 
the particle positions when they are considered 
as stellar (i.e. without collisions; grey) and as 
gaseous (with collisions; black). For this com- 
parison we used / = 0.8, which leads to a mean 
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energy dissipation per particle of 1.25 xlO 3 , 
in good agreem ent with the numbers given by 
ISchward dl98lh . Other numerical values give 
qualitatively similar results. Fig. [9] shows that 
there is in general good agreement between 
the loci of the gaseous and of the stellar arm. 
Furthermore, the gaseous arm is more concen- 
trated than the stellar arm, more so for a larger 
number of collisions per revolution. 

Physically the above results can be under- 
stood as follows. Sticky particles (i.e. gaseous 
particles) follow the same orbits as the stellar 
ones, except at the time of the collisions. This 
ensures a general similarity. Due to the colli- 
sions, however, the gaseous particles lose part 
of their kinetic energy and their velocity ap- 
proaches that of the mean. In Paper I, we com- 
pared the loci of manifolds for different ener- 
gies and found that the ones with the smaller 
energies lie in configuration space within the 
ones with the higher energies. This means that 
the corresponding spiral arms, or rings will be 
thinner for the lower energies. Thus, when the 
particles lose energy they will fall onto an or- 
bit nearer to the mean and the arms will be- 
come thinner. This is exactly what is found 
with the calculations leading to Fig. [9] Thus 
one can think of the lowest energy manifolds 
(i.e. the ones having the energy of L\ and Lj) 
as an attractor, to which the gaseous trajec- 
tories will tend to because of the dissipation 
due to the collisions. Thus the gaseous arms 
should be thinner that the stellar ones, and this 
is indeed observed in real galaxies. This is dis- 
cussed more extensively in Papers III and IV, 
where explicit comparisons are made. 

7. Final remarks 

In this short review and, particularly, in 
the five papers on which it is based 



dRomero-Gomez et all 120061 120071 120081: 



lAthanassoula et alJ 20081 and Athanassoula 



et al. 2009) we presented building blocks 
that can explain the formation of rings (both 
inner and outer) and spirals in barred galaxies. 
The invariant manifolds associated with the 
Lagrangian points L\ and L2 guide chaotic 
orbits which have the right shape for inner 
and outer rings and spirals and thus can be 




Fig. 9. Comparison of the outer ring loci as cal- 
culated with (black) and without (grey) collisions. 
In this example / = 0.8. The positions of the 
Lagrangian points are marked by asterisks. See text 
for a description of the calculations. 



their building blocks. In particular, here we 
discussed morphologies of rR\ type and of 
spiral type. We also discussed the case when 
the L\ and L2 are stable, and four other un- 
stable points bifurcate from them. The whole 
structure then has nine Lagrangian points, five 
stable and four unstable. We also introduced 
an ersatz for gas and showed that the gaseous 
arms and the stellar arms should have roughly 
the same shape and that the former should be 
thinner than the latter. 

Even though dynamically sound and very 
appealing, this theory will be useful for barred 
galaxies only if the manifolds and the orbits 
they guide create structures whose properties 
are in good agreement with observations. This 
point will be fully addressed in paper IV. Here 
let us just underline a couple of noteworthy 
points. 

The first is that the existence of the man- 
ifolds, even if they have the right shape and 
properties, does not necessarily imply that 
the corresponding structures will be present 
in the galaxy. Indeed, the manifolds are only 
the building blocks, and the structure will be 
present only if some mass elements (e.g. stars) 
follow them. This is the same as for periodic 
orbits, which need to trap material around them 
to form the corresponding structure. Whether 
such material will exist or not depends on 
the formation and evolution of the galaxy (see 
Paper IV). 
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The second is that this is not the only the- 
ory attempting to explain spirals and rings, al- 
though most theories attempt either the one or 
the other. Our theory relies on the existence of 
a bar; in other cases a companion is necessary. 
The fact that a theory is correct and reproduces 
well many of the main observational data does 
not necessarily imply that all other theories are 
wrong, or irrelevant. Nor is it necessary in or- 
der to establish a theory to show that all oth- 
ers are wrong. Indeed spirals in different galax- 
ies may have different origins, and even in the 
same galaxy more than one theory can be at 
work. 
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